Nuclear Superfluidity and Cooling Time of Neutron-Star Crust 
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We analyse the effect of neutron superfluidity on the cooling time of inner crust matter in neutron 
stars, in the case of a rapid cooling of the core. The specific heat of the inner crust, which deter- 
mines the thermal response of the crust, is calculated in the framework of HFB approach at finite 
temperature. The calculations are performed with two paring forces chosen to simulate the pairing 
properties of uniform neutron matter corresponding respectively to Gogny-BCS approximation and 
to many-body techniques including polarisation effects. Using a simple model for the heat transport 
across the inner crust, it is shown that the two pairing forces give very different values for the cooling 
time. 



I. INTRODUCTION 

A ncwly-formcd neutron star cools within minutes 
from a temperature of the order of 30 MeV to less than 1 
MeV via neutrino emission. After this stage, the thermal 
evolution of the neutron stars can be strongly influenced 
by the onset of nuclear superfluidity [l|, • This is espe- 
cially the case for rapid cooling models. In these models, 
due to direct Urea or other exotic processes, the core 
cools down so rapidly that a temperature inversion de- 
velops between the core and the crust. The crust acts 
as an insulating blanket which keeps the surface rela- 
tively warm until the cooling wave reaches the surface. 
When this happens, the surface temperature drops pre- 
cipitously to the temperature of the core. One of the 
relevant quantity in this cooling scenario is the cooling 
time, i.e., the time necessary for the cooling wave to ar- 
rive from the cold core to the star's surface. The cooling 
time is primarely determined by the thermal response of 
the inner crust, formed by nuclear clusters immersed in 
a sea of unbound neutrons and ultrarelativistic electrons 

In the rapid cooling models, both the core cooling and 
the time needed for the core-crust thermalisation depend 
critically on nuclear superfluidity. Thus, on one hand, 
the onset of superfluidity in the core matter suppresses 
the neutrino cooling since the total energy of particles in- 
volved in the neutrino production must exceed the pair- 
ing gap. On the other hand, the superfluidity of inner 
crust matter is shortening significantly the cooling time. 
This happens due to the suppression of the heat capac- 
ity of the inner crust matter by the energy gap in the 
excitation spectrum of the superfluid neutron gas. 

One of the first estimation of the cooling time was given 
by Brown et al who considered the possibility of a 
rapid cooling induced by the strangeness condensation. 
They calculated the heat diffusion time through the crust 

where R r is the 



with a simple formula, i.e., tdiff 
thickness of the crust, while C v and n are the specific heat 
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and the thermal conductivity of the crust matter. The 
estimated cooling time was of the order of a few tens of 
years. Later on, using a direct Urea process as cooling 
mechanism, more realistic calculations of thermal evolu- 
tion of neutron stars were performed 0]. The numerical 
simulations showed that the cooling time does not de- 
pend on the details of the rapid cooling mechanism but 
rather on the structure of the neutron star. Besides, it 
was also shown that the cooling time can be strongly re- 
duced (by about a factor of three) if the neutron gas in 
the inner crust is in a superfluid phase. 

In the calculations mentioned above, the effects of nu- 
clear clusters on the superfluid and thermal properties of 
the neutron gas were disregarded. Since then, a few quan- 
tum calculations of the inner crust matter superfluidity, 
including the effects of the nuclear clusters, have been 
done [1,1a, 0|- Thus, using the Hartree-Fock-Bogoliubov 
(HFB) approach it was found that the presence of the nu- 
clear clusters can modify significantly the heat capacity 
of the neutron gas. How the nuclear clusters could affect 
the cooling time of the inner crust was investigated by 
Pizzochero et al § . Using a cooling model similar to the 
one employed by Brown et al. [1| , it was concluded that 
the presence of the clusters, primarily in the outermost 
layers of the inner crust, could change the cooling time 
by amounts comparable with the cooling time itself. It 
was also found that the effect of the clusters on the cool- 
ing time depends rather strongly on the temperature and 
the pairing force used in calculating the specific heat of 
the inner crust matter. 

The impact which the pairing force could have on 
the superfluid properties and the specific heat of the in- 
ner cust matter was recently analysed in the framework 
of HFB approach at finite temperature (FT-HFB) 0]. 
Thus, it was shown that if the pairing force used in the 
FT-HFB equations is adjusted to describe two different 
scenarios for the neutron matter superfluidity, i.e., one 
corresponding to BCS calculations with the Gogny force 
and the other to Gorkov type calculations which take 
into account self-energy and screening effects [9(, the re- 
sults for the specific heat of the inner crust matter can 
change by several orders of magnitude. The scope of the 
present paper is to show what are the consequences of 
these changes in the specific heat upon the cooling time 
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of inner crust matter. In the first part of the paper we 
shall extend the calculations of Ref. Q to the low-density 
region of the inner crust, which was not treated before 
in the HFB approach, and analyse how the specific heat 
and the thermal diffusivity behave across the inner crust. 
Then, using the model of Refs. [4, 8] for the heat trans- 
port, we shall discuss how the cooling time of the inner 
crust depends on neutron matter superfluidity. 

II. THERMAL PROPERTIES OF INNER 
CRUST MATTER IN THE HFB APPROACH 

The thermal response of the inner crust matter de- 
pends on thermal diffusivity, defined as the ratio of the 
thermal conductivity to the heat capacity. The heat ca- 
pacity of the inner crust has contributions from the elec- 
trons, the neutrons and the lattice. The heat capacity 
of the electrons, considered as a uniform and relativistic 
gas, has the standard form [l(| while the contribution of 
the lattice to the specific heat is usually neglected. 

In the normal phase, the specific heat of the neutrons 
exceed the specific heat of the electrons by about two 
orders of magnitude (see Figure 1 below). However, the 
onset of the neutron superfluidity reduces drastically the 
neutron specific heat, which can become smaller than the 
electron specific heat in some regions of the inner crust. 
How the neutron specific heat is affected by the super- 
fluidity as well as by the temperature and the presence 
of nuclear clusters was already studied in Ref. 0], but 
only for a few density regions of the inner crust. Here 
we extend this study to all relevant densities of the inner 
crust, starting from the neutron drip density up to about 
half the nuclear saturation density. This region of the 
inner crust is supposed to give the largest contribution 
to the cooling time of the crust @|. 
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TABLE I: The Wigner-Seitz cells considered in the paper. 
The structure of the cells, i.e., the baryonic densities (p), the 
number of neutrons (N), the number of protons (Z) and the 
cell radii (Rws) correspond to Ref. [llj|. x% are the thickness 
of the layers employed in Eq.(10). 

In microscopic calculations the inner crust matter is 
divided in independent cells treated in Wigner-Seitz ap- 



proximation Up to baryonic densities of the order 
of half the nuclear saturation density, considered in this 
paper, each cell is supposed to contain in its center a 
spherical neutron-rich nucleus surrounded by unbound 
neutrons and immersed in a relativistic electron gas uni- 
formly distributed inside the cell. The proton- to- neu- 
tron ratio and the dimension of the cell at a given bary- 
onic density are determined from the beta equilibrium 
conditions. In the present study we use the cell structure 
determined in Ref. [ll[ by HF type calculations. The 
properties of the cells considered in this paper are dis- 
played in Table I. Compared to Ref. [llj, here we have 
not included the cell with Z=32, which most probably 
belongs to the deformed pasta phase. For the cells listed 
in Table I we shall determine the specific heat by using 
the quasiparticle spectrum generated by the FT-HFB ap- 
proach presented below. 

A. The HFB approach at finite temperature 

The FT-HFB approach for the inner crust matter was 
presented in details in Ref. Q ■ For the sake of complete- 
ness, here we recall the main steps. 

Assuming spherical symmetry for the Wigner-Seitz 
cell, the radial FT-HFB equations have the form: 

h T (r) - A A T (r) \ / U t {r) \ =E ( U^r) \ 
A T (r) -h T {r)+\ ) \ Vi(r) J 1 \ V t (r) J ' 

where Ei is the quasiparticle energy, A is the chemical po- 
tential, hr(r) is the thermal averaged mean field hamil- 
tonian and Ay(r) is the thermal averaged pairing field. 
The latter depends on the average pairing density kt 
given by: 

K T(r) = ^ + m(r)V t (r)(l - 2f t ) , (2) 

i 

where /i = [l+exp(Ei/kBT)]~ 1 is the Fermi distribution, 
ks is the Boltzmann constant and T is the temperature. 
In a self-consistent calculation based on a Skyrme-type 
force, as used here, hxir) depends on the thermal aver- 
aged particle density 

i 

(3) 

as well as on thermal averaged kinetic energy density and 
spin density. The expressions of the last two densities are 
given in Ref. @- 

In the calculations presented here the mean field hamil- 
tonian is calculated with a Skyrme type force while for 
the thermal averaged pairing field we use a density de- 
pendent contact force of the following form [3] : 

V(r-r') = Voll-vi—TWr-r') = V eff (p(r))6(v-r'), 
Po 

(4) 
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where p(r) is the baryonic density and po — 0.16 frn -3 . 
With this force the thermal averaged pairing field is local 
and given by: 



Mr) 



Veff(p(r)) 



K T (r), 



(5) 



where kt{t) is the thermal averaged pairing density. 

To generate in the outer region of the Wigner-Seitz 
cell a constant density corresponding to the neutron gas, 
the FT-HFB equations are solved by imposing Dirichlet- 
Neumann boundary conditions at the edge of the cell 
[Tlj . i.e., all wave functions of even parity vanish and the 
derivatives of odd-parity wave functions vanish. Apart 
from that, the self-consistent solutions of the HF-HFB 
equations are found in the same manner as for finite nu- 
clei. 

The calculation scheme outlined above is employed to 
study how the specific heat of the neutrons is behaving 
in various regions of the inner crust. In order to do that, 
one has to choose the two-body interactions in the FT- 
HFB calculations. These interactions should provide a 
reasonable description of both the nuclear clusters and 
the neutron gas, which are the baryonic components of 
the inner crust matter. For the calculation of the mean 
field we shall use the Skyrme force SLy4 [l3| , which was 
fixed to describe properly the mean field properties of 
neutron-rich nuclei and infinite neutron matter. 

The choice of the pairing force is more problematic 
since at present it is not yet clear what is the strength 
of pairing correlations in neutron matter. Thus, on one 
hand, the BCS calculations with bare forces giv e a maxi- 
mum gap in neutron matter of about 3 MeV [l2| . A max- 
imum gap of about 3 MeV one gets also with the Gogny 
force [151 ] , which is commonly used to describe the pairing 
properties in finite nuclei. On the other hand, if one goes 
beyond the BCS approximation and takes into account 
the in-medium effects, the maximum gap is suppressed. 
The suppression depends on the many-body approxima- 
tions used in the calculations [T^]. In order to analyse 
how the uncertainty on the pairing gap in neutron mat- 
ter could reflect upon the thermal response of the inner 
crust, we shall do calculations with two zero range pair- 
ing interactions which simulate the pairing gap in nuclear 
matter obtained either with the Gogny force, or with 
models which take into account the in-medium effects. 
For the latter we consider a maximum gap of 1 MeV, as 
indicated by recent calculations @. In Ref. [7] the re- 
quirements mentioned above were approximatively satis- 
fied by using two zero range pairing forces (Eq.4) having 
the same parameters for the density dependent term, i.e., 
77=0.7, a=0.45, and two different strengths, i.e., Vq— {- 
430.0,-330.0} MeV fm~ 3 . These values of the strengths 
were obtained by solving the FT-HFB equations with a 
cut-off energy equal to 60 MeV. Since with a 60 MeV 
cut-off we have numerical problems in solving the FT- 
HFB equations for large Wigner-Seitz cells, here we shall 
keep this cut-off and the corresponding strengths only for 
the first two cells while for the other cells we shall take 



a smaller cut-off, equal to 20 MeV. This cut-off is intro- 
duced smoothly, i.e., by an exponential factor e _B *' 100 
acting for quasiparticle energies Ei > 20 MeV. With this 
smooth energy cut-off we shall use the strengths values 
y ={-570.0,-430.0} MeV fm" 3 . The pairing force corre- 
sponding to the first (second) value of the strength will 
be called below the strong (weak) pairing force. 



B. Specific heat 

The quasiparticle spectrum determined by solving the 
FT-HFB equations is used to calculate the specific heat 
of the neutrons inside the Wigner-Seitz cell, i.e., 



_TdS 
Cv ~ Vdf' 



(6) 



where V is the volume of the Wigner-Seitz cell and S is 
the entropy: 

S = -k B ^(2j, + !)(/, In fi + (1 - fi) Ml - /,)). (7) 



The results obtained for the cells listed in Table I are 
shown in Figure 1. In the same figure is also shown the 
specific heat of the electrons, given by [Icl |: 
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,2/3 



3hc 



2/3 



(8) 



The specific heats are calculated for a temperature of 
T=0.1 MeV, which is a typical temperature for the in- 
ner crust matter at the cooling stage analysed here (see 
the discussion below). From Figure 1 we can see that 
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FIG. 1: Specific heat of the neutrons for the Wigner-Seitz 
cells listed in Table I. The results correspond to the strong 
and the weak pairing forces (see the text) and for the cells 
with (without) the nuclear clusters. The specific heat of the 
non-uniform cells obtained when the pairing correlations are 
switched off are indicated by star symbols. The square sym- 
bols show the specific heat of the electrons. 

if the neutrons are in the normal phase, their specific 
heat is greater than the specific heat of the electrons in 
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all Wigner-Seitz cells. When the neutron superfluidity 
is turned on, the specific heat of the neutrons is sup- 
pressed due to the pairing gap in the excitation spectrum. 
Since the suppression depends exponentially on the pair- 
ing gap, the results obtained with the strong and the 
weak pairing forces are very different, as seen for the WS 
cells 1-5. For the second WS cell, in which the pairing 
gap in the neutron gas region have the maximum value, 
the specific heat obtained with the two pairing forces dif- 
fers by about 7 orders of magnitude. In the WS cells 7-10 
the neutron gas is in the normal phase at the tempera- 
ture T=100 keV. Therefore both pairing forces give the 
same results for the specific heat. 

In Figure 1 are shown also the values of the specific 
heat obtained when the nuclear clusters are disregarded. 
For obtaining these values we have just removed the pro- 
tons from the cells and perform the FT-HFB calculations 
in the same conditions as for the cluster+neutron gas. It 
can be seen that in some cases (see the results for the 
cells 2-4) the nuclear clusters could have a sizable influ- 
ence upon the specific heat. However, the influence of 
the nuclear clusters are relatively small compared to the 
effect coming form the uncertainty of the pairing force. 



C. Thermal diffusivity 

The specific heat enters in the heat transport through 
the thermal diffusivity, defined by D — where k is 
the thermal conductivity. In the inner crust, the latter is 
primarly determined by the electrons. The dependence 
of thermal conductivity on density and temperature was 
parametrized by Lattimer et al@, based on the calcula- 
tions of Itoh et al [la]. For a temperature above 10 8 K 
analysed here, the conductivity is nearly independent of 
the temperature and is given by k = C(p/ po) 2 ^ 3 , where 
C = 10 21 ergs cm -1 s _1 . With the conductivity given 



cific heat, the diffusivity is much smaller for the weak 
pairing force, except the last four WS cells. For both 
pairing forces one can see that the diffusivity is much 
smaller in the outermost layers of the inner crust. As 
seen below, these layers have an important contribution 
to the cooling time of the inner crust. 



III. COOLING TIME OF THE INNER CRUST 
MATTER 

In order to calculate the cooling time, i.e., the time 
needed for the cooling wave to propagate from the cold 
core to the surface, one should integrate the heat equa- 
tion 



1 d 2 dT 
r z Or Or 



dT 

~dt' 



(9) 



Since the specific heat and the conductivity depend on 
density and temperature profile of the crust, the solution 
of the heat equation is not trivial. Here we use a sim- 
ple model employed in Refs. 0,11]. The model is based 
on the following assumptions: a) the spherical geometry 
for the heat transport is approximated by a planar ge- 
ometry, i.e, one considers the heat diffusion through a 
one-dimensional piece of matter. This approximation is 
supported by the small thickness of the inner crust com- 
pared to the size of the core; b) the inner crust is divided 
in layers of constant thermal diffusivity. The diffusion 
time through a layer of thickness Xi and diffusivity £>,; is 

calculated by the relation U — l^jj- [HI , where the fac- 
tor 7, which depends on the boundary conditions of the 
problem, is taken equal to 4/pi 2 |8j; c) the total diffu- 
sion time across the crust is obtained by summing up the 
contributions of the layers, i.e., 



D, 



(10) 
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Log( p/p„) 

FIG. 2: Thermal diffusivity (neutrons plus electrons) corre- 
sponding to the Wigner-Seitz cells listed in Table I. The no- 
tations are the same as in Figure 1. 

by this expression and the specific heat calculated in the 
FT-HFB approach one gets the thermal diffusivity shown 
in Figure 2. As expected from the behaviour of the spe- 



In the equation above the thermal diffusivity depends 
on density and temperature, Di = D(p(Ri),T(Ri)) 1 
where Ri is the position of the layer i. In the calculations 
we divide the inner crust into 10 layers, corresponding to 
the 10 cells listed in Table I. The position correspond- 
ing to each cell can be found by solving the Tolman- 
Oppenheimer-Volkov (TOV) equations, which provides 
the density profile of the star. In the present calcula- 
tions we use the solution of TOV equations corresponding 
to the following equations of state [l7|: Baym-Pethick- 
Sutherland [IH for the outer crust, Negele-Vautherin [ll| 
for the inner crust and Glendenning-Moszkowski [l9j for 
the core. From the solution of the TOV equations one 
extracts the radii Ri corresponding to the densities of the 
cells given in Table I. Then, doing a linear interpolation, 
we determine the size Xi of the layers considered for each 
cell. The results are shown in Table I. 

The diffusivity depends also on the temperature pro- 
file. Numerical simulations indicate that before the core- 
crust thermalisation the temperature is increasing from 
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about T=0.1 MeV to about T=0.2-0.3 MeV when one 
goes from the outer part to the inner part of the crust. 
Since the inner part zones of the inner crust have large 
diffusivities, they contribute less to the cooling time com- 
ared to the outermost zones. Therefore, following Ref. 
we shall consider for all layers a flat temperature pro- 
file equal to T=0.1 MeV. The diffusion time across the 
inner crust obtained for this value of the temperature 
is shown in Figure 3. The most striking thing we can 
notice is the critical dependence of the cooling time on 
the pairing force. Thus, for a strong pairing force the 
cooling time is about 12 years. The largest contributions 
come from the outermost zones, as noticed also in Ref. 
Q. Concerning the effect of the clusters, one can see 
that is rather small for this temperature. In the case of 
the weak pairing force, the cooling time is increasing by 
about a factor two compared to the strong force. More- 
over, if the neutron superfluidity is ignored completely, 
the cooling time is further increasing to about 90 years. 
These dramatic changes shows how important is the pre- 
cise knowledge of the neutron matter superfluidity for the 
cooling time of the innner crust. 
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FIG. 3: The diffusion time across the inner crust. The nota- 
tions are the same as in Figure 1. 

A similar strong dependence of the cooling time on the 
pairing scenarios we have obtained by using two other 
zero range forces with the parameters fixed following a 
different protocol, i.e., a unique strength, V"o=-648 MeV 
fm~ 3 , and two sets of parameters for the term depen- 
dend on density, ^={0.95, 0.87 } and a={0.45, 0.2}. The 
value of Vo was taken so that to get, for a smooth cut- 
off energy equal to 20 MeV, the experimental value for 
the scattering length of two free neutrons. The cooling 
times obtained with these two zero range pairing forces 



are equal to about 9.1 and 33.8 years, respectively. 

The cooling times calculated in this section are based 
on the assumption of a flat temperature across the inner 
crust. This is a rather drastic approximation, especially 
for the scenario of a weak pairing force when, as seen in 
Figure 3, all regions of the inner crust contribute signifi- 
cantly to the total diffusion time. More realistic calcula- 
tions of the cooling time should be based on dynamical 
solutions of the heat equations (9). 



IV. SUMMARY AND CONCLUSIONS 

We have estimated the cooling time of the inner crust 
matter using specific heats calculated in the framework 
of HFB approach at finite temperature. In order to study 
the effects of the neutron superfluidity on thermal prop- 
erties of the inner crust, we have employed two paring 
forces. They have been fixed to reproduce the pairing 
properties of infinite neutron matter given either by a 
Gogny force or by microscopic calculations which take 
into acocunt polarisation effects. For the latter we con- 
sidered a maximum pairing gap in neutron matter equal 
to 1 MeV. With the two pairing forces we have studied 
what are the effects of neutron superfluidity on the spe- 
cific heat and the heat difusion of inner crust matter. 
It is shown that the heat difusion predicted by the two 
pairing forces are rather different, especially in the higher 
density part of the inner crust. These differences in the 
heat diffusion have a big influence upon the cooling time. 
Thus, if one shifts from one pairing force to the other the 
cooling time is changing by a factor of three. This show 
how large could be the window in which the cooling time 
may vary due to the present lack of knowleadge of neu- 
tron matter superfluidity. 

The neutron superfluidity affects the cooling time 
through the specific heat, calculated here with the non- 
collective quasiparticle spectrum provided by the FT- 
HFB equations. However, the excitation spectrum of the 
inner crust baryonic matter presents also low-lying col- 
lective modes |20j. Since these modes give an important 
contribution to the specific heat [2l|, they may also af- 
fect significantly the cooling time of the inner crust. This 
issue will be addressed in a future study. 
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